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Using a sample of 88.8 x 10 6 BB events collected with the BABAR detector at the PEP-II storage 
rings at the Stanford Linear Accelerator Center, we measure the branching fractions of seven color- 
suppressed B-meson decays: B(B° -> D°7r°) = (2.9±0.2(stat.) ±0.3(syst.)) x 10~ 4 , B(B° -> ZTV 1 ) 
= (2.9±0.4(stat.)±0.5(syst.)) x 10~ 4 , B(B° -> D°r)) = (2.5±0.2(stat.) ±0.3(syst.)) x 10~ 4 , B(B° -> 
D*°r,) = (2.6±0.4(stat.)±0.4(syst.)) x 10~ 4 , B(B° D°uj) = (3.0 ± 0.3(stat.) ± 0.4(syst.)) x 10~ 4 , 
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B{B° -> D*°uj) = (4.2 ± 0.7(stat.) ± 0.9(syst.)) x 10 -4 , and £>V) = (1.7 ± 0.4(stat.) ± 

0.2(syst.))_x 10~ 4 . We set the 90% confidence-level upper limit: B{B° -> D rjf) < 2.6 x 1(T 4 . The 
channels B° — ► D*°ij, D*°lj, and D ?/ are seen with more than five-sigma statistical significance. 
All of these branching fractions are significantly larger than theoretical expectations based on the 
"naive" factorization model. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

Weak decays like B° -> D^+h~ can proceed through 
the emission of a virtual W~ , which then can materialize 
as a charged hadron [1]. Because the W~ carries no 
color, no exchange of gluons with the rest of the final 
state is required. Such decays are called color-allowed, 
though color- favored might be more apt. By contrast, 
decays like B° — ► D^°h° cannot occur in this fashion. 
The quark from the decay of the virtual W~ must be 
combined with some anti-quark other than its partner 
from the W~. However, other anti-quarks will have the 
right color to make a color singlet only one-third of the 
time. As a result, these decays are "color-suppressed". 
The tree level diagrams for the color-allowed and color- 
suppressed decays are shown in Fig. 1. 

The decays of B° into D^°t: have been observed by 
the CLEO collaboration [2], while the B° decays into 
£)(*)°7T , D°i], D°co, and D°p° have also been measured 
by the Belle collaboration [3, 4]. We present in Table I 
the prior measurements of branching fractions of the B° 
color-allowed and color-suppressed decays. The level of 
color suppression can be estimated from the branching 
fractions for the D^tt and D^p decay modes. 

Since QCD calculations of decay rates from first prin- 
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FIG. 1: The (a) color-allowed and (b) color-suppressed spec- 
tator tree diagrams for B° — > D h decays. 
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TABLE I: Prior measurements of branching fractions for B° 
color-allowed and color-suppressed decays. When two uncer- 
tainties are given, the first uncertainty is statistical and the 
second systematic. We also quote the 90% confidence upper 
limits (UL) when the statistical significance of the measure- 
ment is less than four standard deviations. 



5° mode ff(xnr 4 ) UL (x!0~ 4 ) 



D + ir~ 26.8 ±1.2 ±2.7 [5] 

D°7T 2.9 ±0.5 [6] 

D*+n- 27.6 ±2.1 [6] 

D*°7r° 2.5 ±0.7 [6] 

D + p- 78 ± 14 [6] 

D°p° 2.9 ± 1.0 ± 0.4 [4] 

D* + p- 73 ±15 [6] 

D*°p° - < 5.1 [4] 

D°n l-4±°-l ± 0.3 [3] 

D*°r) 2.0±° l ± 0.4 [3] < 2.6 [6] 

D°oj 1.8± 0.5t u ; 4 [3] 

D*°lu 3.1±{ l ±0.8 [3] < 7.4 [6] 

L>V - < 9.4 [6] 

zry - < 14 [6] 



ciples are at present not possible, we must rely on models 
to describe the above processes. In an early model [7, 8], 
the "naive" (or "generalized" ) factorization model, which 
is very successful in describing charmed meson decays, 
the decay amplitudes of exclusive two-body non-leptonic 
weak decays of heavy flavor mesons are estimated by re- 
placing hadronic matrix elements of four-quark opera- 
tors in the effective weak Hamiltonian by products of 
current matrix elements. These current matrix elements 
are determined in terms of form factors describing the 
transition of the B meson into the meson containing the 
spectator quark, and a factor proportional to a decay con- 
stant describing the creation of a single meson from the 
remaining quark-anti-quark pair. In this approach, the 
decay amplitudes corresponding to Figs. 1(a) and 1(b) 
are proportional to ai and <22 [9], respectively, where the 
di are effective QCD Wilson coefficients. As an example, 
using the naive factorization model, the decay amplitude 
for the B° — > D + n~ mode corresponding to Fig. 1(a) can 
be written as [9] 

A f (B° -> D+tt-) = 

&V cb V: d {ml ml) ai UF B - D (ml), (1) 
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while the decay amplitude for the B° — ► D 7r° mode 
corresponding to Fig. 1(b) can be expressed as [10, 11] 

V2A f (B° ^ D°ir a ) = 

i^V cb V: d (m% - ml)a 2 f D F B ^(m 2 D ), (2) 

where Gp is the Fermi coupling constant, V cb and V ud 
are CKM matrix elements, f v and fo are the decay con- 
stants of the 7T and D mesons, and F$ ~* M {q 2 ) are the 
longitudinal form factors of the B-meson decays to M 
mesons at momentum transfer q 2 . The coefficients a\ 
and ai are real in the absence of final-state interactions 
(FSI) and ideally would be process independent [7-10]. 

The color-allowed B° -> £>(*)+ (tt", p", of) and B -> 

decays, the color-suppressed B — > (cc){K W, 7r) 
decays, and the mixed B~ — > Z?(*)°(7r~, af ) decays 
can all be accommodated by universal constants a\ = 
1.1 ± 0.2 and a 2 ~ 0.2-0.3 [6, 9, 12, 13]. This no longer 
holds for color-suppressed B decays with one c-quark 
only, like D^tt, where measurements listed in Table I are 
inconsistent with a universal value of ai in the absence 
of FSI [10]. The naive factorization model [9, 10, 13-17] 
predicts too small values for the branching fractions of 
the color-suppressed modes, in the range (0.3-1.7) x 10~ 4 
and corresponding to a factor (a2/ai) 2 ~ 0.03-0.09. 

Final state interactions, however, may change this pic- 
ture significantly and, thus, may increase substantially 
these rates, as rescattering effects can connect the final 
states shown in Fig. 1(a) and Fig. 1(b) (see, for example, 
Ref. [16]). Similar effects have already in the past com- 
pletely changed the conclusions of the models that de- 
scribe non-leptonic D° decays, especially for decay modes 
such as D° -> K°ir° [18]. Therefore, in the case of large 
FSI, a description in terms of isospin amplitudes is more 
appropriate and will be used in Sec. IX B to discuss our 
results. 

This situation is an impetus for higher precision re- 
sults and the investigation of additional channels that 
might provide clues to the underlying mechanisms. In 
this paper we report on the branching fraction measure- 
ments of the seven color-suppressed B°-meson decays to 
£>(*)°7r°, D(*)°2, DWV and D°rf . We also report on a 
search for the B° — > D*°rf decay. These results are based 
upon an integrated luminosity equivalent to 88.8 x 10 6 
BB events. This corresponds to about nine times that 
used for the earlier measurement by CLEO [2] (9.7 x 10 6 
BB events) and about four times that used for the ear- 
lier measurements by Belle [3] (23.1 x 10 6 BB events). 
Recently, with 31.3 x 10 6 BB events, the Belle collab- 
oration has reported branching fraction measurements 
for B° — > Z?(*)°7r + 7r~ decays, including the D°p° mode, 
as already discussed, and the investigation of the D*°p° 
channel [4]. We present the first measurement of the 
B° — > D*°r), D*°ui, and D°r]' modes with more than 
five-sigma statistical significance. 



II. THE BABAR DETECTOR AND DATA 
SAMPLE 



The BABAR detector is located at the PEP-II e+e~ 
storage rings operating at the Stanford Linear Acceler- 
ator Center. At PEP-II 9.0-GeV electrons collide with 
3.1-GeV positrons to produce a center-of-mass energy of 
10.58 GeV, the mass of the Y(4S). The data used in this 
analysis were collected with the BABAR detector and cor- 
respond to an integrated luminosity of 81.9 fb _1 recorded 
at the T(45) resonance. 

The BABAR detector is described in detail in Ref. [19]. 
Surrounding the interaction point is a 5-layer double- 
sided silicon vertex tracker (SVT), which gives precision 
spatial information in three dimensions for charged parti- 
cles and measures their energy loss (dE/dx). The SVT is 
the primary detection device for low-momentum charged 
particles. Outside the SVT, a 40-layer drift chamber 
(DCH) provides measurements of the polar angles and of 
the transverse momentum (pt) of charged particles with 
respect to the beam direction, together with the SVT. 
The resolution of the px measurement for tracks with mo- 
menta above 1 GeV/c is up T /pT = 0.13% x p T + 0.45%, 
where px is measured in GeV/c. The drift chamber mea- 
sures dE/dx with a precision of 7.5%. Beyond the outer 
radius of the DCH is a detector of internally reflected 
Chercnkov radiation (DIRC), which is used primarily for 
charged-hadron identification. The detector consists of 
quartz bars in which Cherenkov light is produced when 
relativistic charged particles traverse the material. The 
light is internally reflected along the length of the bar 
into a water-filled volume mounted on one end of the 
detector. The Cherenkov rings expand in the water vol- 
ume and are measured with an array of photomultiplicr 
tubes mounted on its outer surface. A CsI(Tl) crys- 
tal electromagnetic calorimeter (EMC) is used to detect 
photons and neutral hadrons, as well as to identify elec- 
trons. The resolution of the calorimeter can be expressed 
as a E /E = 2.3%/(E)i 1.9%, where E is measured in 
GeV. The EMC detects photons with energies down to 
20MeV. The EMC is surrounded by a superconducting 
solenoid, which produces at 1.5-T magnetic field. The in- 
strumented flux-return (IFR) consists of multiple layers 
of resistive plate chambers (RPC) interleaved with the 
flux-return iron. The IFR is used in the identification of 
muons and long-lived neutral hadrons. 

Signal and generic background Monte Carlo events 
are generated using the BABAR particle decay simulation 
package [20], the "EvtGen" package. The interactions of 
the generated particles traversing the detector are sim- 
ulated using the GEANT4 [21] program. Beam-induced 
backgrounds, which varied from one data-taking period 
to the next, are taken into account in the simulation of 
the detector response. This is done by adding the sig- 
nals generated by these beam-induced backgrounds to 
the simulation of the various physics events. 
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III. PARTICLE RECONSTRUCTION AND 
COUNTING OF BB EVENTS 

Charged-particle tracks are reconstructed from mea- 
surements in the SVT and/or the DCH. The tracks 
must have at least 12 hits in the DCH and px > 
100 MeV/c [22]. In the case of the tracks used to recon- 
struct mesons, we also use tracks reconstructed with 
the SVT alone (see Sec. IV B 1). The tracks must extrap- 
olate to within 20 mm of the e + e~ interaction point in the 
plane transverse to the beam axis and to within 50 mm 
along the beam axis. Charged- kaon candidates are iden- 
tified using a likelihood function that combines dE/dx 
and DIRC information. The likelihood function is used 
to define tight and loose kaon criteria as pion vetos. To 
satisfy the tight kaon criterion, the track must also have 
p > 250 MeV/c and make an angle with respect to the 
electron beam direction, which is used as the reference 
axis for all the polar angles, between 0.45 and 2.50 rad 
so that the candidate is within the fiducial region of the 
DIRC. Photons are identified by energy deposits in con- 
tiguous crystals in the EMC. Each photon must have an 
energy greater than 30 MeV and a lateral shower shape 
consistent with that of an electromagnetic shower. 

The measurement of branching fractions depends upon 
an accurate measurement of the number of BB meson 
pairs in the data sample. We find the number of BB pairs 
by comparing the rate of spherical multi-hadron events 
in data recorded on the Y(4S) resonance to that in data 
taken off-resonance. This latter data sample is collected 
40 MeV below the T(4S) resonance and corresponds to 
an integrated luminosity of about 10 fb" 1 . 

The purity of the multi-hadrons events is enhanced by 
requiring the events to pass selection criteria based on all 
tracks (including those reconstructed in the SVT only), 
detected in the fiducial region 0.41 < 6 < 2.54 rad and 
on neutral clusters with an energy greater than 30 MeV, 
in the fiducial region 0.410 < 6 < 2.409 rad: 

• There must be at least three tracks in the fiducial 
region. The total energy of the charged and neutral 
particles in the fiducial region must be greater than 
4.5 GeV. 

• The ratio of the second to the zeroth Fox- Wolfram 
moment [23] must be less than 0.5. All tracks and 
neutral clusters defined above are used. 

• The event vertex must be within 5 mm of the nom- 
inal beam-spot position in the plane transverse to 
the beam and within 60 mm along the beam direc- 
tion. 

These requirements are about 95.4% efficient for BB 
events as estimated from Monte Carlo simulation. The 
systematic uncertainty on the number of BB events is 
1.1%. 



IV. MESON CANDIDATE SELECTION 
A. General considerations 

The color-suppressed B° meson decay modes are re- 
constructed from D Q or D* Q meson candidates that are 
combined with light neutral-meson candidates h° (jr , T], 
u), and r)'). Events are required to pass the selection cri- 
teria used for BB counting listed in Sec. III. Additional 
requirements discussed below are applied to the signal 
sample. 

We combine tracks and/or neutral clusters to form can- 
didates for the mesons produced in the B decays. Vertex 
constraints are applied to charged daughters before com- 
puting their invariant masses. At each step in the decay 
chain we require that mesons have masses consistent with 
their assumed particle type. If daughter particles are 
produced in the decay of a parent meson with a natural 
width that is small relative to the reconstructed width, 
we constrain the meson's mass to its nominal value. This 
fitting technique improves the resolution of the energy 
and the momentum of the B° candidates as they are 
calculated from improved energies and momenta of the 
DW° and h°. 

We select D*°, D°, h°, and B° candidates using only 
well-understood discriminating variables in order to re- 
duce the systematic uncertainties for the branching frac- 
tion measurements. We choose selection criteria that 
maximize the quality factor Q = S/ + B, where S 
and B are the expected number of signal and background 
events. The values of S and B are estimated from signal 
and background Monte Carlo simulation and data in the 
signal sidebands, but not from data in the signal regions. 
When optimizing the cuts, the values of S have been es- 
timated using the previous branching fraction measure- 
ments obtained by the CLEO [2] and Belle [3] collabo- 
rations. For the D^°rj' analyses, a conservative value 
for the branching fractions equal to 10~ 4 has been as- 
sumed. In most cases we find that Q does not change 
significantly when selection criteria are varied near their 
optimal values. This allows us to choose selection criteria 
that are common to most final states. 

B. Selection of h° and p ± candidates 

The momentum of the h° candidate must satisfy the 
condition 1.3 < p* < 3.0 GeV/c. This requirement is 
loose enough that various sources of background populate 
the sidebands of the signal region. These sidebands are 
used in the background estimate for the signal. 

1. 7T° and p ± selection 

The 7T° meson is reconstructed from photon pairs. We 
consider three sources of ir° with decreasing momenta: ir° 
originating from B° decays, from D°, 77, and u> decays, 
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and directly from D*° decays. The latter two sources are 
discussed below. The mass resolution of ir° candidates 
from B° decays with momenta p* near 2 GcV/c is dom- 
inated by the uncertainty in the opening angle between 
the two photons and is approximately 8 MeV/c 2 . 

These 7r°s are also combined with charged pions to at- 
tempt the reconstruction of p~ mesons. The charged 
pions arc not required to satisfy our regular selection 
criteria for tracks. Thus we retain also low momen- 
tum charged pions that are reconstructed with the SVT 
alone. A 7r 7r~ pair is selected if its mass is recon- 
structed within 250 MeV/c 2 of the nominal p~ meson 
mass. The p~ candidates are used to reconstruct the 
color-allowed B~ — > D^°p~ decays that form a signif- 
icant background for B° — > D^°tt°. The color-allowed 
decays have branching fractions about fifty times that 
for B° -> DM°7r° and the ;y mimic the latter through an 
asymmetric p~ decay in which the tt° carries most of the 
available energy. We veto events with a reconstructed 
B~ — > D(*~)°p~. A discussion of the veto is deferred until 
Sees. VIA and B. 



2. r\ selection 

The r\ candidate is reconstructed in the 77 and 
7r + 7T~7r° decay modes. The branching fraction in the 77 
mode is almost twice as large as that of the 7r + 7r~7T° de- 
cay channel and the efficiency for the 77 mode is greater 
since there are fewer particles to detect. 

In the 77 decay mode we require that the photons have 
energies greater than 200 MeV. A photon is not used if 
it can be paired with another photon with energy greater 
than 150 MeV to form a ir° candidate with an invariant 
mass in the range 120-150 MeV/c 2 . The mass resolution 
for rj — > 77 is approximately 15 MeV/c . 

In the 7r + 7T~7r decay mode, the 77 meson is recon- 
structed employing a vertex constraint that requires a 
X 2 probability greater than 0.1%. To reduce combinato- 
rial background the charged-pion candidates must have 
momentum greater than 250 MeV/c and they must fail 
the tight kaon criterion, while the 7r° must have an en- 
ergy greater than 300 MeV and a mass in the range 1 15— 
150 MeV/c 2 . The mass resolution for rj — > 7r + 7r~7r° is 
approximately 4 MeV/c 2 . 

3. u> selection 

The u> meson is reconstructed in its ir + ir~ir decay 
mode, employing a vertex constraint that requires a \ 2 
probability greater than 0.1%. To reduce combinatorial 
background, the charged pion candidates must have mo- 
mentum greater than 200 MeV/c and they must fail the 
tight kaon criterion, while the ir° must have an energy 
greater than 250 MeV and a mass in the range 120- 
150 MeV/c 2 . The mass resolution of the u> is dominated 
by its natural width of approximately 10MeV/c 2 . The 



use of additional angular properties in the u> meson de- 
cays will be described in Sec. IV D 1. 



4- rf' selection 

We reconstruct the rf meson in its ir + ir~ r/(-^ 77) de- 
cay mode. The product of the branching fractions of sec- 
ondary decays in this channel is 17.5% [6] . This limits the 
signal efficiency, so a separate event selection for D^°rf 
is used. We use the n + TT^rj decay mode rather than the 
dominant p ^ mode as it provides a much cleaner signal. 

The two photons used to reconstruct the r\ candidate 
are required to have energies greater than 100 MeV. A 
photon is not used to reconstruct the r\ meson if it can 
be paired with another photon with energy greater than 
100 MeV to form a ir° candidate with mass in the range 
120-150 MeV/c 2 . We select r\ candidates with a mass 
in the range 495-600 MeV/c 2 . To obtain the highest 
possible signal efficiency we rely on the high purity of the 
signal and impose neither a momentum nor any particle- 
identification requirement on the charged pions. For the 
same reason, a vertex constraint is applied to the 7r + 7r~ 
pair when computing the energy and the momentum of 
an rf meson candidate, but there is no requirement on 
the x 2 probability of the vertex. The mass resolution for 
rf — ► TT + -K~ri(-^> 77) is approximately 4 MeV/c 2 . 



C. Selection of D° and D*° candidates 

The momentum of the Z)(*)° mesons must satisfy the 
condition p* > 1.5 GeV/c. As for the light neutral- 
hadron selection, this requirement retains sidebands, 
which can be used to evaluate backgrounds. 



1. D° -> K~n + , K-tt+tt , and Tr+Tr+Tr" 
selection 

The D° mesons are reconstructed in three decay 
modes: K~n + , K~ir + ir°, and K~ir + ir + n~ . The x 2 
probability for the vertex fit of the charged pions is re- 
quired to be greater than 0.1%. In the K~tt + final state 
the kaon candidate must satisfy the pion veto require- 
ment, while in the K~n + ir° and K~Tr + ir + Tr~ final states 
the kaon candidate must satisfy the tight kaon criterion 
because of the increased background present in these 
combinations. All pion candidates must fail the tight 
kaon criterion. 

To reduce combinatorial background in the K~n + Tr a 
final state we use the results of the Fermilab E691 exper- 
iment [24], which determined the distribution of events 
in the Dalitz plot. This distribution is dominated by the 
two possible K* resonances (K* a — ► K~ir + or K*~ — > 
K~ir°) and by the p + {-^ 7r + 7r°) resonance. We select 
only those events that fall in the enhanced regions of the 
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Dalitz plot as determined by experiment E691. Recon- 
structed 7T° mesons are required to have masses in the 
range 115-150 MeV/c 2 . The mass resolution is approx- 
imately 6.5 MeV/c 2 . To increase the signal purity only 
7T° mesons with energy greater than 300 MeV, as defined 
in the laboratory frame, are retained. 

The D° mass resolutions are approximately 6.7, 
10.7, and 5.0 MeV/c 2 for the K~-k+, K-jt+tt , and 
K~ir + iT + Tr~ decay modes, respectively. 



2. D*° -f D°n° selection 

The D*° mesons are reconstructed in the D°ir° de- 
cay mode. The D° candidates are selected as described 
above. The ir° candidates are required to have momenta 
that satisfy the condition 70 < p* < 300 MeV/c and a 
mass in the range 115-150 MeV/c 2 . The mass resolution 
for the soft ir° daughter is approximately 6.5 MeV/c 2 . 
The resolution of the D*°~D° mass difference is approx- 
imately 1 MeV/c 2 . 



D. Selection of B candidates 

1. Event shape and angular distributions 

Both BB events and u, d, s, and c quark- antiquark 
events contribute to the combinatorial background that 
does not peak near the nominal B mass. To reject u, d, 
s, and c components we use shape variables and angular 
distributions that distinguish these from the signal BB 
events. 

Because the u, d, s, and c continuum events are jet-like, 
while B meson decays produce spherical events, we can 
suppress them by requiring that the ratio of the second to 
the zeroth Fox- Wolfram moment [23] must be less than 
0.5 as described in Sec. III. For each reconstructed B° 
candidate we compute the thrust and sphericity axes of 
both the candidate and the rest of the event, using only 
the tracks and neutral clusters as defined in Sec. III. We 
define the angles 6 t hr and 6 sp h between the axes of the B° 
candidate and the rest of the event. The distributions of 
|cos6>thr| and |cos# sp h| peak near 1.0 for u, d, s, and c 
background while they are nearly flat for B decays. Thus 
we require at least one of the conditions | cos# sp h| < 0.85 
or |cos6» thr | < 0.85 to be true for the D^tt , D^°r], 
and Z?M°w modes. Since the two angles $thr and $ sp h are 
strongly but not completely correlated for signal events, 
the relative signal efficiency for this requirement is close 
to 92%. This is larger than the relative signal efficiency 
of about 85% if only the requirement |cos# t hr| < 0.85 
is applied, while the background rejection is about the 
same. 

For the L>(*)°7r°, D^°r], and D^°lj final states we 
also take advantage of the sin 2 6b* distribution of the 
polar angle 6b* ■ This quantity is the angle between the 
B momentum vector and the beam axis in the T(4S*) 



rest frame. Therefore we only keep the candidates that 
satisfy | cos 6b* \ < 0.8 as the distribution is almost flat 
in | cos 6b* | for combinatorial background. 

For the D^°rf channels, we have seen that the event 
yield is expected to be small. In order to keep the sig- 
nal acceptance as high as possible, we use a more com- 
plex scheme. We require |cos# t hr| < 0.9 and then cal- 
culate a Fisher discriminant (T) that combines eleven 
variables [25]. Two of these are the two polar angles 6b* 
and 9t, where 9t is the angle between the B candidate 
thrust axis and the beam axis in the T(4S) rest frame. 
The other nine are the scalar sums of the energies of all 
charged tracks and neutral showers (except those used 
in the B candidate reconstruction) binned in nine 10° 
polar angle intervals relative to the B candidate thrust 
axis. The separation between the means of the signal and 
qq background distributions of the T variable is 1.2—1.3 
times the width of either distribution. 

For the D°u> channel where the uj is necessarily longi- 
tudinally polarized, we use the properties of the distribu- 
tions of two additional angles. The angle #n is the angle 
between the normal to the plane of the three daughter 
pions in the u> center-of-mass frame and the line-of-flight 
of the i3-meson in the u rest frame. The angle #d is 
the angle, in the rest frame of one dipion, between the 
third pion and either of the other two. The signal events 
are distributed as cos 2 6>n and sin 2 On , while the corre- 
sponding cos6>n and cos6*d distributions are nearly flat 
for combinatorial background. We select only events in 
a region of the three-dimensional parameter space of the 
angles 0b* , #n , and #d that has high signal efficiency. 
This region is defined by 



cos6» D 2 + cos6> B . 2 < 0.64, 



and 



cos 6*D 
0.8 



cos 6b* 
0.8 



+ 



cos6» N | - 1.0 
05 



cos#n| — 1-0 
05 



< 1, 



< 1. 



(3) 



(4) 



(5) 



In the D*°lo channel, the u polarization is not known 
a priori and we apply only the requirement given by 
Eq. (3). 

For the D*°h°, h° = ir°, i], and i]' modes where the 
D*° is longitudinally polarized, we use the angular decay 
distribution to reject combinatorial background. The an- 
gle 6>hd is defined as the angle between the line-of-flight 
of the D° and the one of the Z?*°, both evaluated in the 
D*° rest frame. The distribution is almost flat in cos 6>hd 
for combinatorial background, while signal events are dis- 
tributed as cos 2 #hei- Fot the D*°tt° and D*°r] channels 
we require 



cos 6b 
0.8 



+ 



cosflhcil -1.0 
0.6 



< 1. 



(6) 
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For the D*°f]' final state we only require | cos#hci| > 0.4 
since the angle 6*b* is already included in the definition 
of T. 



2. Multiple B candidates 

After applying the above selection criteria, a small 
fraction of events have more than one B candidate. The 
average multiplicity of B candidates for the data events is 
between 1.01 and 1.18, depending on the D° decay mode. 
The average multiplicity is slightly higher for the D*°h° 
modes than for the D°h° modes. With the exception of 
the D^°r]' final states we select the B candidate with 
the lowest value of 
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where a mD and a mh are the resolutions of the measured 
D° and h° masses. The last term in the equation is only 
present for D*° decays and a/\ mB , D is the average resolu- 
tion of the measured D*°-D° mass difference. The mass 
resolutions depend on the decay modes and are slightly 
different for data and Monte Carlo simulation. Each of 
the three terms is found to be approximately Gaussian 
with mean value near zero and standard deviation near 
one. 

In order to reduce combinatorial backgrounds, we re- 
quire that each of the terms in Eq. (7) is less than 2.5 2 
(a ± 2.5ct requirement). In the case of the u mesons, 
the candidates must have a reconstructed invariant mass 
within 25 MeV/c 2 (± 2.5 times the uj natural width) of 
the nominal value. 

For the £)(*)°jj' channels, the signal acceptance is rel- 
atively lower than for other modes, but the background 
level is also much smaller. Therefore we keep all the can- 
didates in the events and weight them by 1/N where N 
is the number of B candidates in the event. In order to 
reduce the combinatorial background for these two chan- 
nels the invariant mass of the rj' candidate is required to 
be within 2.5 a of its nominal value. The D° candidates 
are required to have a reconstructed mass within 2-3 a 
(depending on the decay mode) of their nominal value. 
We reject D*° candidates whose D*°-D° mass difference 
is not within 3 a of its nominal value. 



3. B candidates and background yields 

Two kinematic variables arc used to isolate the B- 
meson signal for all modes. One is toes, the beam- 
energy-substituted mass. The other is AE, the differ- 
ence between the reconstructed energy of the B candi- 
date and the beam energy in the e + e~ center-of-mass 



frame. Both quantities use the strong constraint given 
by the precisely known beam energy (the beam energy 
is known to within a fraction of an MeV). The beam- 
energy-substituted mass is defined as 



' s/2 + pp.pB 
En 



(8) 



and the energy difference is 

AE = E* D + E* h - y/l/2. (9) 

Where y/s is the e + e~ ccntcr-of-mass energy. The small 
variations of the beam energy over the duration of the 
run are taken into account when calculating toes- For 
the momentum pi (i = 0,B) and the energy Eq, the 
subscripts and B refer to the e+e~ system and the 
reconstructed B meson, respectively. The energies E^ 
and El are calculated from the measured £)(*)° and h° 
momenta. Signal events have toes — m B a and AE ~ 0, 
within their respective resolutions. 

We limit the selection of the B° candidates to the 
"signal neighborhood", defined by |A1?|< 350 MeV and 
5.2 < to es < 5.3 GeV/c 2 . The toes resolution is dom- 
inated by the beam energy spread and is approximately 
3MeV/c 2 , depending slightly on the B decay mode. 
The AE resolution for the D^°ir° and D^°t](-^ 77) 
modes is dominated by the angular and energy resolu- 
tion of the EMC. The AE resolution is approximately 
37-44 MeV for the L>(*)°7r° modes and 28-35 MeV for 
the £)(*)°?j(— ► 77) modes, depending on the D*° and 
D° decay mode. The AE resolution is better for the 
D°r](-> tt+tt-tt ), £>W°w, and D«V modes because 
the angular and the momentum resolution for charged 
tracks is better than for photons. For these modes it is 
approximately 15-20 MeV. 

We define the signal region using the resolutions in 
Toes and AE obtained from the Monte Carlo. The limits 
of the signal region are 5.270 < toes < 5.290 GeV/c 2 
(about ± 3 a around the B mass) and \AE\ < 3 a. 
In the case of B° — * D^ n decay modes, we reduce 
the contribution from the color-allowed B~ — ► D^°p~ 
background by requiring AE to be in the region from 
—90 to 100 MeV. We change these requirements slightly 
for the D^°r]' channels where we want to optimize the 
statistical significance. Here the signal region is defined 
by \AE\ < 2-3 a depending on the D° decay mode and 
5.273 < toes < 5.286 GeV/c 2 . The number of signal 
candidates is computed in the signal region for each B° 
decay mode and the signal Monte Carlo simulation is 
used to determine the acceptance. 

We perform an unbinned maximum likelihood (ML) fit 
to the toes distribution to extract the number of signal 
candidates (-A/" C and)- A fit to the toes distribution allows 
us to model the signal and background shapes with a 
well known, simple, and universal function, independent 
of the B decay mode analysed. 

In the fit the signal component is modeled by a Gaus- 
sian distribution whose a is constrained to the value ob- 
tained from the signal Monte Carlo separately for each 
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TABLE II: The number of candidates (A/ an d), the number of non-peaking (A/" np b) and peaking (A/" p b) background events, the 
number of cross-feed (A/cf) background events from other color-suppressed modes, the number of signal events (S) after peaking 
and cross- feed backgrounds are subtracted, and the statistical significance of the signals (S/^JS + Abckgd)- We obtain A/" np b 
from a fit to the data ?71es distribution, while A/pb is estimated from the Monte Carlo simulation. The statistical uncertainty 
on S includes the uncertainty on N can d as obtained from the ML m E s fit. The statistical uncertainty on A/pb and the estimated 
uncertainties for Acf are accounted for in the systematic uncertainties of the branching fractions. For the Z)'*'°»/ modes, the 
number of can didates is small; t herefore Poisson statistics rather than Gaussian statistics are used. The statistical significance 
is defined as ^/ '2\n{C m ax / 'jC(O)) , where C m ax is the likelihood at the nominal signal yield and £(0) is the likelihood with the 
signal yield set to 0. In the table, the symbol "-" means that the corresponding number can be neglected. 



B° mode 
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A/"npb 


A/"pb 


A^CF 




5 


statistical 


(decay channel) 














significance 


D°TT° 


556 ± 34 


603 ± 22 
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487 ± 34 


14.3 


D*°7r° 


102 ± 12 


32 ± 6 


11 ± 5 


2 ± 


1 


88 ± 12 


7.6 


D°V(-^ 77) 


200 ± 20 


181 ± 12 


17 ± 3 


10 ± 


2 


173 ± 20 


8.9 




76 ± 12 
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1 


74 ± 12 


6.2 


D*°V(-^ 77) 


43 ± 7 


8 ± 2 




4 ± 


1 


40 ± 7 


5.5 


D°uj 


207 ± 18 


136 ± 10 


4 ± 3 


5 ± 


1 


198 ± 18 


10.7 


D* a u> 
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58 ± 7 




5 ± 


1 


70 ± 12 


6.1 


D° V ' 


27 ± 6 


10 ± 1 








27 ± 6 


6.3 




4 ± 2 










4 ± 2 
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B° decay mode. The value of A/ ca nd is computed from the 
fit within the m E s signal region defined earlier. The back- 
ground component is modeled by an empirical phase- 
space distribution [26] (henceforth referred to as the AR- 
GUS distribution): 

A(mES]mo,£,,a) = a tuesV^ — (wes/iwo) 2 x 

exp(£(l - (m ES /mo) 2 )), (10) 

where Too is set to a typical beam energy (5.29 GeV), a 
is the fitted normalization parameter, and £ is the fitted 
parameter describing the shape of the function. 

The ML fit is performed within the limits of the signal 
region in AE, as defined above, and for toes between 
5.2 and 5.3GcV/c 2 . For the D^°r}' modes, in addition 
to using the Toes resolution obtained from the Monte 
Carlo simulation, the mean of the Gaussian distribution 
is also constrained in the ML fit to the nominal B mass. 
The value of the £ parameter in the ARGUS function is 
fixed to the value obtained from a ML fit to the ibes 
data in the AE sideband: 200 < \AE\ < 350 MeV and 
5.2 < m ES < 5.3GeV/c 2 . 

The ARGUS function accounts for random combina- 
torial background originating from u, d, s, and c contin- 
uum events, t + t~ events, two-photon processes, and BB 
events but not for "peaking background" from B°B° and 
B + B~ decays, which have distributions that peak in the 
same location as signal events do. The number of non- 
peaking-background events (A/ np b) is determined from 
the fit to the data in the full 5.2 < m ES < 5.3 GeV/c 2 
interval and the AE signal region by integrating the AR- 
GUS function over the much smaller signal region. 

The number of peaking-background events (A/pb) is 



small relative to the non-peaking background but it is 
dangerous because the peaking-background events lie in 
the signal region. Peaking background comes also from 
color-suppressed decays in B°B° events that are incor- 
rectly reconstructed. This small contribution (A/cf) is 
evaluated separately and thus does not contribute to the 
value of A/pb, as discussed in Sec. V. Altogether we write 
the total number of background events (A/,kgd) in the 
signal region as 

A/bkgd = A/npb + A/pb + A/cf • (11) 
Finally, the number of signal events is calculated as 

S = A/and _ A/pb - A/cF ■ (12) 

The values of A/and, A/" n pb, A/" p b, A/cf, S, and the statis- 
tical significance of the signals for the B° decay channels 
studied in this paper are listed in Table II. 

V. BACKGROUND ESTIMATION 

A. Peaking backgrounds from BB decays other 
than color-suppressed modes 

To investigate backgrounds that peak at the B mass in 
the toes distribution, we use two types of Monte Carlo 
samples: a sample that contains only B~ — > £)(*)°p _ 
and a sample that contains all other charged and neutral 
B-meson decays, except the color-suppressed B° decay 
modes reported in this paper. In the next section we 
describe how we estimate the cross-feed from the color- 
suppressed B° decay modes. 
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The peaking background is estimated with a ML fit 
to the Monte Carlo samples, using a Gaussian distribu- 
tion for signal and an ARGUS background distribution, 
just as for the data (see Sec. IV D 3). We constrain the 
ARGUS shape parameter £ to be the same as the one ob- 
tained for the corresponding data toes distribution. The 
normalization of the ARGUS function is a free parameter 
as are all parameters of the Gaussian. The values of the 
parameters of the Gaussian distribution for the peaking- 
background events are expected to be different than that 
for signal events. The mean value of the Gaussian dis- 
tribution is possibly different from the B mass and the 
resolution is expected to be larger than the nominal value 
for signal events, which is about 3MeV/c 2 . 

The peaking background is taken to be the area under 
the Gaussian distribution in the signal region 5.270 < 
toes < 5.290 GeV/c 2 (5.273 < to es < 5.286 GeV/c 2 
for £)(*)°iy' channels), normalized to the luminosity of 
the data. Table II gives the estimate of the number of 
peaking-background events to be subtracted from the fit- 
ted candidate event yields in the data for each of the 
various channels. For each channel, the number is the 
sum of the various contributions estimated from the BB 
background Monte Carlo samples. As this number is ex- 
tracted from Monte Carlo simulations, we use the statis- 
tical uncertainty associated with this quantity as a sys- 
tematic uncertainty for the branching fraction measure- 
ments. 

The systematic uncertainty due to the constraint ap- 
plied to the ARGUS parameter £, which is fixed to the 
data value in the ML fit to the various Monte Carlo Toes 
distributions used for the peaking-background computa- 
tion, is small or negligible. This systematic uncertainty is 
estimated by recalculating the peaking background when 
using two other fixed values for £. These two values 
are computed from ML fits to two mss distributions ob- 
tained with the Monte Carlo simulation. One distribu- 
tion corresponds to the sum of all the normalized con- 
tributions from the various background sources (peaking 
or non-peaking) only. The second one also includes the 
expected contribution from the signal events. It is found 
that the values of £ for the two types of Monte Carlo 
Toes distributions are very close (within the statistical 
uncertainties) to the corresponding data value. 



B. Peaking backgrounds from other 
color-suppressed modes 

Signal event yields must be corrected for cross-feed be- 
tween color-suppressed modes. Cross-feed occurs when a 
true decay chain of type k is erroneously reconstructed 
as a candidate decay chain of type j. This will bias the 
signal yield for events of type j if such events of type k 
enter the signal region. Cross-feed to each signal from 
B° — > D^°h° decays is investigated using signal Monte 
Carlo samples for these decay modes. In the end, we find 
that the contribution of cross-feed is for the most part 



less than half the statistical uncertainty in the signal. 

For each light neutral hadron type, ft , the dominant 
contribution to B° —> D°(D*°)h° arises from the associ- 
ated B" -» D*°(D°)h° mode. In the case of the D*°h° 
decay modes, since we only consider the D*° — ► D°7r° 
channel, the contribution from the final state D*° —> 
is non negligible. These cross-feed contributions peak at 
the same toes as the signal, but are shifted in AE. 

The number N k ^j of events of type k entering the 
signal region for type j is given by 



N k ^j = N(BB)B k A k . 



(13) 



where N(BB) is the number of BB pairs and B k is the 
branching fraction of the decay chain k including the B° 
branching fraction. A k ^j denotes the probability for an 
event of type k to enter the signal region for decay mode 
j. The probability A k ^j is estimated from the Monte 
Carlo simulation as 



A k — yj 



SlAC.k- 



gen,fc 



(14) 



Here, SMC.fc^j is the number of events of type k entering 
the signal region for decay mode j and N gen ^ k is the num- 
ber of generated Monte Carlo events. It is convenient to 
introduce the fractional cross-feed quantity 
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N k ^ 3 _ B k Ak- 



BjAj^j 



(15) 



For a given candidate event of type j, the probability 
that it is generated by one of the possible cross-feed con- 
tributions can be expressed by the fraction FcfU) given 
by 



or, using Eq. 15, by 



1 + E w ^ 



(16) 



(17) 



In what follows, FcfU) is simply written as Tcf for each 
color-suppressed decay mode j. 

In order to calculate N k ^j , we must know the branch- 
ing fractions of the investigated decay modes. We use 
recently measured values for the branching fractions of 
the ft , D°, and D*° decays chains [6]. We consider 
63 color-suppressed B° — ► D^°h° decay chains. The 
light neutral hadron ft is a n°(^> 77), an > 77 
or 7T + 7r~7r°), an ui(— > 7T + 7r~7r°), a p°(-^ 7r + 7r - ), or an 
»?'(—» n + -K~ri(-^> 77) or p°-f) meson. The D° mesons 
are reconstructed in the modes K~tt + , K~n + Tr a , and 
K ~tt + tt + it~ , and the D*° mesons in the channels D°7r° 
and D°7. For the B° — > D^°h° branching fractions 
we use the values measured in this analysis (summarized 
in Table VIII). These final branching fractions are de- 
termined after several iterations because the cross-feed 
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TABLE III: Total fractional cross-feed {T) expressed in per- 
cent (see text for definition) observed in the Monte Carlo 
simulation. The dominant sources that contribute are shown 
in decreasing order of importance. 



B° mode J : cf(%) Dominant sources 



D°ir° 


3.6 


D*°TV° 


D*\D%°)^ 
7A(77) 


2.6 


D*\D°j)7v , D°tv° 
D*°ri, D°tt°, D*°tt° 


5.4 


DVtt+tt-tt ) 


2.2 


D*° V , D°tv°, D°lo 


L»* (D 7r°)r;(77) 


8.8 


D*°(D° 7 ) V , D*°7T , D° V 


D°LJ 


2.5 


D*°uj 


D* (D°ir°)uj 


6.5 


Zr (D° 7 )a;, D°cu, 
and D°r;(7r + 7r-7r ) 







estimate depends upon the branching fractions being 
measured. Therefore, we iterate the calculation of the 
background from cross-feed until the values of the com- 
puted branching fractions do not change by more than 
1CT 6 . For the contributions from D°p° and D*°p° chan- 
nels we use the results obtained recently by Belle [4]: 
(2.9 ± 1.0(stat.) ±0.4(syst.)) x 1CT 4 and the upper limit 
5.1 x 10~ 4 , respectively. In the latter case the assump- 
tion of such a large value for the branching fraction is 
likely to be an overestimate; yet the Z)(*)°p° decays do 
not generate any significant cross-feed contributions to 
any of the modes studied in this paper. 

Table III shows the total contributions from cross-feed 
to each mode reported in this study. The dominant 
sources are also shown in decreasing order of importance. 
The number of cross-feed events, A/qf, is calculated as 
the difference between the number of candidates in the 
data and the number of other peaking-background events 
estimated from the Monte Carlo simulation, which in- 
cludes no signal, multiplied by the fractional cross-feed: 

■A/cf = (Arcane! - W p b) x T C f- (18) 

The corresponding number of cross-feed events is listed 
in Table II for each mode. 

The cross-feed contributions for the B° — ► D^°rf 
analyses are found to be negligible. This is due to both 
the good mass resolution of the mode r\' — > ir + ir~i](-^> 
77) and to the complexity of the signature used to re- 
construct these signals. 

VI. B° CANDIDATES IN THE VARIOUS 
COLOR-SUPPRESSED DECAY MODES 

A. B° D°ir° 

Figures 2(a) and 2(b) show the distributions in toes 
with -90 < AE < 100 MeV and in_A£ with 5.270 < 
toes < 5.290 GeV/c 2 for candidate B° -> D°n° events. 



The solid line in Fig. 2(a) represents the ML fit to the 
sum of the ARGUS and Gaussian functions. In Fig. 2(b) 
the hatched histograms represent the simulated events 
for the signal and separately for the various backgrounds 
from BB and qq (q = u, d, s, c) events. 

Peaking backgrounds originate from color-allowed 
B~ — > D°p~ decays where the n~ from the p~ — ► -K~n° 
decay has very low momentum and is missed in the re- 
construction of the D°7T final state. This type of back- 
ground populates the AE plot in the region that is at 
least one pion mass below the signal region. It produces 
a peak in the toes distribution in and slightly below the 
signal region. Resolution effects in AE will cause some 
events to migrate from below the signal region into the 
signal region and thus contribute to the signal peak in 
the toes distribution. 

A veto on the color-allowed B~ — > D°p~ decays is 
applied as part of the selection of the B° candidates. A 
B° candidate is rejected if it can be reconstructed as a 
B~ — > D°p~ candidate with the following properties: 

• It uses the same D° and w° as the B° candidate and 
the p~ meson is selected as described in Sec. IV B 1. 

• The toes is within 9 MeV/c 2 of the nominal B~ 
mass and \AE\ < 100 MeV. 

According to the Monte Carlo simulation, this veto re- 
moves only a few percent of signal events, while it rejects 
about 70% of D°p~ events and 60% of D*°p~ events. 
This background reduction occurs nearly entirely in the 
AE region below approximately one pion mass and the 
veto is less effective in the signal region, where only a few 
percents of the background events are rejected. 

The veto is nevertheless very useful because it de- 
creases the AE distribution in the region just below the 
signal region, thereby reducing the likelihood that the fi- 
nite energy resolution will shift events from the negative 
AE region into the signal region. The precise determi- 
nation of the resolution here is related to the resolution 
of the EMC for relatively energetic ir° mesons. Remov- 
ing a large fraction of these background events at and 
below the lower signal region limit reduces substantially 
this uncertainty. Even after the veto is applied, as it can 
be seen in Fig. 2(b), the shape of the AE distribution for 
this background changes abruptly at about minus one 
pion mass and that below this limit the magnitude of the 
B~ -» D^°p~ background can still not be neglected. 

The yield of the fitted candidate D°ir° events and the 
numbers for the various background contributions to this 
decay mode are listed in Table II. 

B. B° -> D*°ir° 

Figures 2(c) and 2(d) show the distributions in toes 
with -90 < AE < 100 MeV and in AE with 5.270 < 
TO ES < 5.290 GeV/c 2 for the candidate B° -> D* n° 
events. 
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FIG. 2: Distributions of toes and of AE for (a, b) candidate B° —> D°tv° events and (c, d) candidate B° — > _D*°7r° events. 
The dots with error bars correspond to data. In the mBS distribution, the ARGUS and Gaussian ML fits are superimposed. 
The number of signal candidates (M can d), which includes peaking-background and cross-feed contributions, is the area of the 
Gaussian function in the signal region 5.270 < m ES < 5.290 GeV/c 2 . The non-peaking background (A/" np b) is represented by 
the shaded region. The hatched histograms in the AE distributions represent the simulated events, and are shown separately 
for signal and the various backgrounds from BB and qq (q — u, d, s, c) events. 



The D*°7r° candidates are contaminated by color- 
allowed B~ — > flW ' p~ decays. Events from B~ — > 
D(*)°p- can enter the signal region when the soft 7r~ 
from the p~ decay is missed. In the case of D°p~ events 
an unrelated 7r° is used to reconstruct the D*° meson. 
For this mode we veto both B~ — > D°p~ and D*°p~ 
decays. The criteria used to veto B° candidates are the 
same as for the veto described in the D°7r° subsection ex- 
cept that for D*°tt° the B° candidate is rejected if there 
is a B~ — > £>(*'°p~ candidate that uses the same D" 
and ir° mesons as the B° candidate. 

According to the Monte Carlo simulation this veto re- 
jects about 65% of D°p~ events and 70% of D*°p~ events 
and the signal efficiency is close to 80%. The veto is rel- 
atively less effective in the signal region of the AE dis- 
tribution, where 60% of D°p~ events and 40% of D*°p~ 



events are rejected. As in the D 7r° case discussed above, 
the veto reduces the systematic uncertainty related to the 
background estimate. 

The yield of the fitted candidate D*°ir° events and the 
numbers for the various background contributions to this 
decay mode are listed in Table II. 



C. B° — D°r, 

Figures 3(a) and 3(b) show the distributions in tues 
with \AE\ < 89MeV (3 times the AE resolution mea- 
sured in the Monte Carlo simulation) and in AE with 
5.270 < toes < 5.290 GeV/c 2 for candidate B° -> D°i] 
events, where the rj meson is reconstructed in the 77 de- 
cay channel. Figures 3(c) and 3(d) show the same distri- 
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FIG. 3: Distributions of tties and of AE for (a, b) candidate B° —+ D°r/ (77 — + 77) events, (c, d) candidate B° — » D°r/ 
(77 — > iy + tv~iy°) events, and (e, f) candidate B° — > D*°rj (r) — > 77) events. The various contributions are shown as in Fig. 2. 



butions when the 77 meson is reconstructed in the 7r + 7r~7r° 
decay channel. Here the selection |AU| < 54MeV is ap- 
plied (again 3 times the AE resolution) in the mss dis- 
tribution. 

In the 77—5-77 case, the contribution to the peaking 
background from B~ — > D^*>°p~ decays is dominant. It 
corresponds to 80% of the peaking background. In this 
case a photon from the fast 7r° in the p~ decay is com- 
bined with another photon to form an 77 candidate. This 
background is sufficiently suppressed by the n° veto de- 
scribed in Sec. IV B 2 so that no additional requirements 
are imposed. 

According to the Monte Carlo simulation, the peaking 
background is negligible in the 7r + 7r~7r° decay channel. 
The Monte Carlo simulation includes processes such as 
D^K-n-TT+n that may fake a B° ^ D°rj{-^ n+ir-n ) 
signal if one charged 7r is lost in the reconstruction of the 
B° meson. The branching fractions for these modes have 
been measured recently by the CLEO collaboration [27]. 
Because these backgrounds are shifted in AE by more 



than the mass of the missing n and because the 77 mass 
selection is quite tight, the Monte Carlo simulation in- 
dicates that no events originating from such modes are 
selected within the signal region. We checked the effect 
of widening the signal region to \AE\ < 110 MeV. Due 
to resolution effects more background events in the AE 
sideband region migrate into the wider AE signal region; 
we observe that in that latter case about 10% of the to- 
tal BB background is generated by DW?7(— » 7r + 7r~7r°)7r - 
decays. 

The yields of the fitted candidate D°t] events for the 
V ~ * 77 an d 7T + 7T~7r decay modes and the numbers 
for the various background contributions to these decay 
modes are listed in Table II. 



D. B° -> D*°T7 

Figures 3(e) and 3(f) show the distributions in mES 
with I AS I < 92 MeV (3 times the AE resolution mea- 



17 




FIG. 4: Distributions of m,ES and of AE for (a, b) candidate B° — > D°uj events and (c, d) candidate B° — > D*°ui events. The 
various contributions are shown as in Fig. 2. 



sured in the Monte Carlo simulation) and in AE with 
5.270 < m ES < 5.290 GeV/c 2 for candidate B° -> D*°7] 
events in which the r\ meson is reconstructed in the 77 
channel. 

According to the Monte Carlo simulation, the peaking 
background is negligible. The yield of the fitted candi- 
date D*°rj events and the numbers for the various back- 
ground contributions to this decay mode are listed in 
Table II. The statistical significance of the signal is 5.5. 

E. B° -> D°u> 

Figures 4(a) and 4(b) show the distributions in toes 
with \AE\ < 61MeV (3 times the AE resolution mea- 
sured in the Monte Carlo simulation), and in AE with 
5.270 < toes < 5.290 GeV/c 2 for the candidate B° -> 
D a u> events. Due to the tight u> mass selection and the 
angular selections, the peaking background is small. 

For the peaking-background determination, we have 



included possible contributions from D^*' it" ir~ tt + de- 
cays. CLEO [27] reports the observation of these pro- 
cesses, gives branching fractions for D*7r~7r~7r + 7r° and 
D(*'u>w~, and provides evidence for D^*'p'~(-^> ujtt^). 
These measurements have been performed for both 
charged and neutral B decays. If the additional 7r _ from 
the p'~ decay is missed, these decays can fake B° — ► D°ui 
events. But Monte Carlo simulation indicates that the 
AE distribution for this background is shifted by more 
than the mass of the missing pion and rarely falls in the 
signal region. We estimate from the Monte Carlo simula- 
tion that about 11% of the total B°B° background in the 
signal region originates from D^ + ujtt^ modes; similarly, 
13% of the total B + B~ background is from DWW- 
decays. These fractions remain the same if the AE sig- 
nal range is extended to \AE\ < 100 MeV, thus indi- 
cating that the DW 7t~ 7t~ tt + n° background is randomly 
distributed in AE over the signal region. We also find 
that D(*'p events contribute about 5% of the total BB 
background. 
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FIG. 5: Distributions of tties and of AE of (a, b) candidate B° — » D°i]' events and (c, d) candidate B° — > D*°rf events. 



The yield of the fitted candidate D°u> events and the 
numbers for the various background contributions to this 
decay mode are listed in Table II. 



F. B° -> D*°u 

Figures 4(c) and 4(d) show the distributions in mES 
with \AE\ < 61MeV (3 times the AE resolution mea- 
sured in the Monte Carlo simulation) and in AE with 
5.270 < m ES < 5.290 GeV/c 2 for candidate B° -> D*°lu 
events. 

As for the D°ui analysis, when determining the peak- 
ing background, the effect of tt~ tt~ tt + tt° decays has 
been evaluated. In this case the mode B~ — > D°cutt~ 
may contaminate the signal when the ir~ is replaced by 
a 7T° to fake a D*° meson. However, the kinematics of 
the soft tt° in the D*° decay for the D*°u> signal is very 
different from those of the tt~ where the momentum can 
be large. In addition the relatively small branching frac- 
tion for the D°u>ir~ decays implies that the contribution 



from this background in the signal region is not expected 
to be important. We estimate from the Monte Carlo sim- 
ulation that 16% of the total B + B~ background in the 
signal region originates from D^°lu7t^ modes. No con- 
tribution to the B°B° background from the D^*^ + o;7r~ 
decays has been found. The fractions remain the same 
when the AE range of the signal region is extended to 
\AE\ < 100 MeV. Again, this confirms that this type 
of background is uniformly distributed in AE over the 
signal region and rules out any significant contribution 
to the peaking background from these decays. We also 
find that D^p events contribute about 5% of the total 
BB background. The peaking background for this decay 
mode is found to be negligible. 



The yield of the fitted candidate D* a ui events and the 
numbers for the various background contributions to this 
decay mode are listed in Table II. The statistical signifi- 
cance of the signal is 6.1. 
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TABLE IV: Acceptance (.4), corrected acceptance (Ax>rr) ob- 
tained after differences between Monte Carlo simulation of 
detector response and data are taken into account, and over- 
all efficiency {£ ) that includes branching fractions from sec- 
ondary decays. The uncertainties associated with these num- 
bers are discussed in Sec. VIII. 



B° mode 


A (%) 


Aorr (%) 


£(%) 


(decay channel) 








D°ir° 


9.1 


7.9 


1.87 


1)*%° 


2.7 


2.3 


0.34 


D° V (-^ 77) 


9.7 


8.6 


0.82 


D° V (-^ 7V + 7Y-TY ) 


6.5 


5.6 


0.30 


D*°V(-^ 77) 


3.3 


2.8 


0.17 


4.2 


3.5 


0.75 


D*°U! 


1.7 


1.4 


0.19 


D%' 


5.0 


4.2 


0.18 


D*°T]' 


1.6 


1.4 


0.035 



G. 



7T+7T and 



Figure 5(a) shows the toes distribution for D° 
K-TT+ with \AE\ < 60MeV and for D° -» K 
K-TT+tr+n- with \AE\ < 40MeV. Figure 5(b) shows 
the AE distribution with 5.273 < m ES < 5.286 GeV/c 2 . 
According to the Monte Carlo simulation, the peaking 
background in this channel is negligible. As reported in 
Table II, the fit yields Af cand = 26.6 ± 6.0 candidate D°r)' 
events and A/" n pb = 10.4 ± 1.1 combinatorial-background 
events. The statistical significance of the signal, calcu- 
lated from Poisson statistics, is 6.3. 



H. 



D° 



and 



Figure 5(c) shows the toes distribution for 
K--K+ with \ AE\ < 60MeV and for D° -> K~ 
K-tt+tt+tt- with \AE\ < 40MeV. Figure 5(d) shows 
the AE distribution with 5.273 < to es < 5.286 GeV/c 2 . 
According to the Monte Carlo simulation the peaking 
background is negligible. As reported in Table II, the 
fit yields Af ca nd = 4.0 ± 2.2 candidate D*°rf' events and 
A/"n P b = 0.5 ± 0.3 combinatorial-background events. The 
statistical significance of the signal, calculated from Pois- 
son statistics, is only 3.0. 



VII. BRANCHING FRACTIONS 

The acceptance A for signal events is estimated from 
signal Monte Carlo as 



A = 



Smc 



N, 



(19) 



TABLE V: Values of Aorr, B(D°) (the branching fraction 
of the various D° decay modes [6]), £> sec (the product of the 
branching fractions associated with the secondary decays of 
the 7? — » 77 and the D°), and £ for the B° — ► D°?7(77) decay 
mode. The branching fraction for the r\ — > 77 is taken to be 
39.4% [6]. The uncertainties associated with these numbers 
are discussed in Sec. VIII. 



D° 


decay 


^4corr (9o) 


B(D°)(%) 


B sec (%) 


£(%) 


K~ 


7V + 


19.5 


3.8 


1.5 


0.29 


K~ 


7V+7V" 


6.0 


13.1 


5.1 


0.31 


K~ 




7.4 


7.5 


2.9 


0.22 


all 




8.6 




9.5 


0.82 



Where Smc is the number of events in the signal region 
that pass the selection criteria and N gen is the number of 
generated signal Monte Carlo events. 

The selection efficiencies for each mode are obtained 
from detailed Monte Carlo studies in which the detector 
response is simulated using the GEANT4 [21] program. 
The efficiencies of tracking, detection and reconstruction 
in the EMC, vertex fitting, and particle identification 
have been measured in control sets of data and compared 
with their Monte Carlo simulation. We correct the accep- 
tance for differences between data and Monte Carlo sim- 
ulation of these effects by using precise correction factors 
that are applied to each track (for track reconstruction 
efficiency), to each photon, 7r°, 77(77) (f° r neutral cluster 
detection efficiency and energy resolution), to each kaon 
candidate (for particle identification efficiency), and to 
each vertex- fit (for vertex- fit efficiency). Most of these 
corrections depend upon the polar angle and momenta 
of the tracks and neutral clusters and some also depend 
on the running conditions. 

Tracking efficiencies are determined by identifying 
tracks in the SVT and measuring the fraction of tracks 
that are reconstructed in the DCH. The 7 and ir° efficien- 
cies are measured by comparing the ratio of the number 
of events N(t+ -» V T h+TT°) and N{t+ -» t7 T /i+7r 7r°) to 
the known branching fractions [28]. The kaon identifi- 
cation efficiency is estimated from a sample of D* + — > 
D°ir+, D° -> K-TT+ decays that are identified kinc- 
matically. Based on a similar selection, a sample of 

B° -» D* + 7T-, D*+ -» D°7T+, D° -» K-7T+, K-JT+TT , 

or if~7r + 7r + 7r~ decays is used to determine the vertex-fit 
efficiency corrections. 

The acceptances A obtained with Eq. (19) and the 
corrected acceptances A CO n are listed in Table IV. The 
last column in Table IV lists the values of the overall 
efficiency £ defined as 
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FIG. 6: Measured branching fractions for each of the three 
D° decay modes and for the combination of the three for (a) 
B° -> D°7T and (b) B° -» D*°7r°. The shaded bands rep- 
resent the results from the present investigation. The length 
of the error bars is equal to the sum in quadrature of the 
statistical and the systematic uncertainty; the statistical con- 
tribution is superimposed on the error bar. The CLEO [2] 
and Belle [3] results are also shown. 



where 



B scc = B(D*° -> £>°7r°) x S(tt° -> 77) x 
B(h° ->Y)x J2 B ( D ° -» x ) 



(21) 



v 



is the product of the branching fractions associated with 
the secondary decays of the D*°, h°, and D° (with 
X = K--K+, K-TT+n", or K-tt+tt+tt-). The B{D*° -> 
D°7r°) x B(tt° — > 77) factor is only present for the 
B° — > D*°h° final states. Note that the overall efficiency 
£ for the D^*'°rf decays is reduced with respect to the 
other B° modes by the relatively small values of B scc . 

In Table V we display, as an example, the contributions 
of the three D° final states in the decay mode B° — > 
_D°?7(77). There are variations between the acceptance 
and branching fraction for the three D° decay modes 
leading to similar values of £ for the three modes. A 
similar conclusion holds for other B° — > D^°h° final 
states. 

To obtain branching fractions, the number of back- 
ground subtracted signal events, S, is divided by the 
number of BB events in the data sample, N(BB), and 
the overall efficiency, £: 



B(B a -> D^°h°) = 



S 



N(BB) x £ 



(22) 



These branching fraction calculations assume equal pro- 
duction of B°B° and B + B~ pairs at the T(AS) reso- 
nance. 



VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are associated with the accep- 
tance corrections discussed in Sec. VII. The uncertain- 
ties from the tracking-efficiency corrections are 0.8% per 
charged track. To take into account uncertainties caused 
by the vertex reconstruction, we assign a systematic un- 
certainty equal to 1.1% per two-track vertex and 2.2% 
per four-track vertex. For particle identification the un- 
certainty is 2.5% per track. The uncertainties from 
the requirement that all the daughters must fail the 
tight kaon criterion are negligible. Uncertainties in the 
acceptances for photon detection account for imperfect 
simulation of photon-energy and position resolution, thus 
affecting ir° and rj reconstruction efficiencies and the AE 
resolution. For the detection of isolated 7r° and 77(77) 
mesons uncertainties of 5% and 2.5% are used. These 
uncertainties are summed in quadrature, together with 
other corrections that depend upon the energy of each 7 
used to reconstruct the mesons. 

We consider systematic uncertainties from other 
sources. For the cross-feed fractions an uncertainty equal 
to 25% of the estimated fraction accounts for uncertain- 
ties in the branching fractions reported in this study and 
used in the cross-feed determination. This value is cho- 
sen conservatively; it corresponds to the branching frac- 
tion measurement with the largest uncertainty reported 
in this paper (see Table VIII). 

The effect of the specific AE range used to define the 
signal region and based on the resolution measured from 
the Monte Carlo simulation has been estimated by vary- 
ing the limits of the range by ± 0.5 a. The observed 
variations in the branching fraction are used to deter- 
mine the systematic uncertainty from this source. In the 
case of the D^^ir modes, we vary the lower limit on the 
signal region definition (—90 < AE < 100 McV) between 
— 110 and —60 MeV. Therefore, this procedure also ac- 
counts for uncertainties in the peaking-background esti- 
mates that are caused by the systematic uncertainty of 
the energy resolution that originates from the EMC. 

To evaluate the systematic uncertainty associated with 
using the wes resolution taken from Monte Carlo in the 
fit to data, we also let it vary freely in that fit and half of 
the variation in the yields is taken as the systematic error. 
We also investigate the uncertainties in the combinatorial 
background due to setting the value of the ARGUS shape 
parameter £ to the value obtained in the fit to the data 
toes distribution in the upper AE sideband + 6 a < 
AE < 350 MeV. For the Z)(*) ?/ analyses, the value of £ 
is obtained from the AE sidebands (see Sec. IV D 3). We 
therefore vary the value of £ by one standard deviation 
of the statistical error. In each case we take half the 
variation observed as the systematic uncertainty. Finally, 
the sum of the systematic errors from the ARGUS shape 
parameter and the fixed Gaussian width is taken as the 
systematic error for the toes fitting procedure. 

Systematic uncertainties in the peaking background 
arise from imprecision in branching fractions and from 
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TABLE VI: Systematic uncertainties of the measured branching fractions in percent, 
systematic uncertainty is negligible. 



The symbol "-" indicates that the 



Category 


D 7T° 


£>* 7T° 


DV77) 


D°?7(7r + 7r-7r ) 


D*°r) 


D°u 


D*°u) 


D°f)' 


D*°ri' 


Tracking 


2.1 


2.1 


2.0 


3.6 


2.0 


3.6 


3.6 


3.6 


3.6 


Vertex-fit 


1.4 


1.4 


1.4 


2.5 


1.4 


2.5 


2.5 


1.4 


1.4 


Kaon identification 


z.o 


z.o 


z.o 


z.o 


z.o 


z.o 


z.o 


O EC 
Z.O 


z.o 


7, 7T°, and r\ detection 


5.2 


8.1 


3.7 


6.0 


6.8 


5.9 


9.1 


3.5 


6.5 


Cross- feed 


l.U 


U. ( 


1 /I 


u.o 


9 /I 


u.o 


1 7 






AE resolution 


1.7 


1.9 


3.0 


4.4 


3.5 


5.7 


3.3 






m-ES fit 


0.3 


3.2 


4.5 


4.8 


8.4 


3.0 


10.3 


2.3 


2.3 


Peaking background 


3.3 


6.3 


3.2 


2.0 


0.5 


3.4 


4.0 






Event selection 


6.8 


9.4 


6.1 


8.9 


7.6 


6.8 


11.9 


7.9 


7.9 


B(D W °) and B(h°) 


4.6 


6.6 


4.4 


4.6 


6.3 


4.3 


6.4 


5.6 


7.3 


Number of BB pairs 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Monte Carlo statistics 


0.7 


2.2 


1.3 


1.6 


2.0 


1.6 


2.9 


1.6 


1.6 


Total (%) 


11.1 


16.4 


11.2 


14.5 


15.8 


13.5 


20.8 


11.7 


13.7 



statistical uncertainties in the number of peaking- 
background events obtained from the procedure de- 
scribed in Sec. VA. For the Z?(*)°7r° modes the system- 
atic uncertainty associated with the veto of the B~ — > 
D^°p~ background has been studied and is part of the 
systematic uncertainty of the background estimate. For 
these £)(*)°7r° decay modes, we remove the veto on the 
B~ — > D^°p~ background and we include in the un- 
certainties half of the relative variation of the branching 
fraction. Finally we have explained in Sec. VA how the 
systematic uncertainty related to the fitting method used 
in the calculation of the number of pcaking-background 
events is estimated. The variation of the branching frac- 
tion due to the latter effect is small or negligible (4% at 
most) but is included in the systematic uncertainty from 
peaking background. 

We vary the selection criteria applied to several other 
uncorrelated variables such as invariant masses, event 
shape, and helicity angles (see Sees. IV B, IV C, and 
IV D). We conservatively assign a single systematic un- 
certainty due to the efficiencies associated with these 
many selection criteria, equal to the quadratic sum of the 
average of the absolute values of the observed changes in 
branching fraction for each variable. None of the various 
observed changes contribute in a dominant way to the 
total systematic uncertainty due to event selection. 

The uncertainties from the counting of BB pairs, from 
the branching fractions of and h° secondary de- 

cays [6] , and from the statistics of the Monte Carlo sam- 
ples used to determine the signal acceptance, are also 
considered. 

The systematic uncertainties described above are listed 
in Table VI for all the modes reported in this paper. It is 
seen that the dominant systematic uncertainties are due 
to the event selection, from 7, 7r°, and r\ detection, from 
the ttjes fitting procedure, and from the £>M and h° 



branching fractions. 



IX. RESULTS 
A. Branching fractions 

The branching fractions of the color-suppressed modes 
reported in this paper and their statistical and system- 
atic uncertainties are listed in Table VII for the three D° 
decay modes K~n + 7 if~7r + 7r°, and K~Tr + ir + TT~ . The 
measurements obtained by combining the three D° de- 
cay modes are presented in Table VIII. Except for the 
D*°r]' decay channel all measurements have statistical 
significance in excess of five-standard deviations. For the 
D*°i]' decay channel we quote a 90% confidence level up- 
per limit using Poisson statistics. To aid in combining our 
result with future results for D*°r]' a central value with 
statistical and systematic uncertainties is also given. For 
the D°rj decay mode the most precise result is obtained 
by combining the r\ — > 77 and r\ — > 7r + 7r~7r° decay modes. 

The results listed in Tables VII and VIII are also pre- 
sented in the summary Figs. 6, 7, and 8 for comparison. 
It is seen that, for a given B decay, the three measure- 
ments using the three D° decay modes are consistent 
among themselves. Where available, previous results by 
the CLEO [2] and Belle [3] experiments are also shown. 
The precision of the results on the branching fractions 
presented in this paper can be compared to the precision 
of existing measurements as listed in Tab. I. 

In some cases theoretical predictions are more precise 
for ratios of branching fractions than for branching frac- 
tions themselves [15, 29]. An example is the ratio of 
B(B° -» £>(*)V) to B(B° -» DM /?) [15]. Systematic 
uncertainties partly cancel in the measurement of ratios 
so they are also more precisely determined experimen- 
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TABLE VII: Measured branching fractions for B° — > D^*'°h° (xlO 4 ). The measurements are given for each of the three D' 
decay modes K~tt + , K~it + tt , and K~tt + tt + it~ . The first uncertainty is statistical and the second systematic. 
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FIG. 7: Measured branching fractions for each of the three 
75° decay modes and for the combination of the three for (a) 
B° -> D°rj(^ 77), (b) B° -> D°rj(^ tt+tt-tt ), (c) each of 
the B° — » D°rj modes and their combination, and (d) B° — > 
D*°r]. The branching fraction for the D°r) channel is obtained 
as the average of the branching fraction of each of the two 
r\ decay modes, weighted by the statistical uncertainties of 
these decays. The computation of the systematic uncertainty 
includes both the correlated and uncorrelated uncertainties of 
these two modes. The Belle [3] results are also shown. The 
error bars are as in Fig. 6. 



FIG. 8: Measured branching fractions for each of the three 
D° decay modes and for the combination of the three for (a) 
B° -> D°u, (b) B° D*°u, and (c) B° -> D°rj' . The 
Belle [3] results, when existing, are also shown. The errors 
bars are as in Fig. 6. 
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TABLE VIII: Measured branching fractions for B° -> D { * )0 h a 
obtained by combining the three D° decay modes. The first 
uncertainty is statistical and the second systematic. The last 
column is the statistical significance. The branching fraction 
for the D°r] mode is obtained as the average of the branch- 
ing fractions of each of the two rj decay modes, weighted by 
the statistical uncertainties of these decays; the computation 
of the systematic uncertainty includes both the correlated 
and uncorrelated errors of these two modes. For the D^°i]' 
modes, the number of candidates is small, so Poisson statistics 
rather than Gaussian statistics are us ed, and the value for the 
statistical significance is defined as y/2\n(£ rnax /£(0)), where 
Cmax is the likelihood at the nominal signal yield and £(0) 
is the likelihood with the signal yield set to 0. For the D*°ij' 
decay mode we also quote a 90% confidence level upper limit 
using Poisson statistics. 



B° mode 


B (xl(T 4 ) 


statistical 


(decay channel) 


significance 


D°ir° 


2.9 ± 0.2 ± 0.3 


> 6.5 


D*°TT° 


2.9 ± 0.4 ± 0.5 


> 6.5 


D° V (-^ 77) 


2.4 ± 0.3 ± 0.3 


> 6.5 


D° V {-^ tv+tv-tv ) 


2.8 ± 0.4 ± 0.4 
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D°ri (combined) 


2.5 ± 0.2 ± 0.3 
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D°U! 
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4.2 ± 0.7 ± 0.9 


6.1 
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1.3 ± 0.7 ± 0.2 
< 2.6 (90% CL) 


3.0 



TABLE IX: Ratios of branching fractions for B° -> D ( * )0 h°. 
The first uncertainty is statistical, the second is systematic. 



B ratio 



this experiment theoretical prediction 



B ^D* uj 



B ^D*°n' 

B ^D°r,' 
B"^D° V 

B°^D*°r)' 



1.0 ± 0.1 ± 0.2 1.20 [16] or 1.0 ± 0.2 [29] 

0.78 [16] 



0.9 ± 0.2 ± 0.1 



0.7 ± 0.1 ± 0.1 



1.3 ± 0.8 ± 0.2 



0.7 ± 0.2 ± 0.1 



0.5 ± 0.3 ± 0.1 



0.41 [9] 
0.64-0.78 [15] 
0.64-0.68 [15] 
0.67-0.68 [15] 



TABLE X: Comparison of our measurements and theoretical 
predictions of naive factorization for the branching fractions 
of the color-suppressed B° decays reported in this paper. For 
the experimental results the statistical and systematic uncer- 
tainties are added in quadrature. 



B° decay mode 


B(xW- 4 ) 


S(xl0- 4 ) 




this experiment 


theory (factorization) 


D°ir° 


2.9 ±0.4 


0.58 [16] (0.7 [9]) 


D*%° 


2.9 ±0.6 


0.65 [16] (1.0 [9]) 


D% 


2.5 ±0.4 


0.34 [16] (0.5 [9]) 


D*°r, 


2.6 ±0.6 


0.37 [16] (0.6 [9]) 


D°uj 


3.0 ±0.5 


0.66 [16] (0.7 [9]) 


D*°oj 


4.2 ± 1.1 


1.7 [9] 


D° V ' 


1.7 ±0.4 


0.30-0.32 [15] 


D*°r)' 


1.3 ±0.7 


0.41-0.47 [15] 



tally. We compare measured ratios of branching ratios 
to theoretical predictions in Table IX. 



B. Isospin symmetry and decay amplitudes 

Isospin symmetry relates the amplitudes for the B~ — > 
D(*)°7r-, B° -> £)(*)+7r-, and B° -» D^V decay 
modes [14]. These amplitudes can be expressed as [10] 



A(D^°7r-) = V3A 3/ ^ DM , 
A(DW+tt-) = y/T/3A 3/2iDW + y/V3A 1/2iDM , 
V2A(D^°7r°) = v / 473^3/2,c(«) - #^1/2,00, 



and 



(23) 



where the amplitudes Ai/ 2 . d(*) and A 3 / 2 ,dW correspond 
to pure 7 = 3/2 and 7 = 1/2 isospin eigenstates. This 
leads to the triangle relation: 

A(D^°n-) = A(D^+n-) ± V2 A(D^°tt°). (24) 

If the relative strong-interaction phase between the two- 
isospin amplitudes {5 D (*)) is equal to zero, the interfer- 
ence between these isospin eigenstates is maximally de- 
structive for the color-suppressed B° — > D^°h° decay, 
while it is respectively maximally constructive for the 
color-allowed B a -> D^+h~ decay. It follows from QCD 
factorization [13], in the heavy-quark limit, that 

/V2A 3/2iDM = 1 + 0(A QC D/m Q ), (25) 

where tuq represents m c or m& and where the correction 
to "1" is also suppressed by a power of 1/N C , the number 
of colors [10, 30]. The above relation also implies that 
S D (,) = C(Aqcd/^q)- Final-state interactions (FSI) ef- 
fects in the 7 = 3/2 and I = 1/2 channels might be 
expected to be independent, leading to a non-zero phase 
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difference 5 D (*> . If the value of 5 D (,) is large enough it 
will substantially undo the destructive interference for 
the color-suppressed decay B° — ► D^°h°, increasing the 
associated branching fraction. 

Using the various equations listed above, the val- 
ues from Table VIII for B(B° -> D^ir ), the Par- 
ticle Data Group values [6] for B(B~ — > D*°ir~) and 
B(B° — > Z?*+7r~), the recent measurements by the CLEO 
collaboration [5] for -» D°n-) and 2?(B 

£> + 7r~), and the B meson lifetime ratio t(B + )/t(B°) = 
1.083 ± 0.017 [6], we calculate the value of the strong 
phase difference \Sd\ = 30° ± 5° for Dir final states and 
\5d* I = 33° ±5° for D*tt final states. The ratio of isospin 
amplitudes \Ai/ 2 .dm / V^A 3 / 2 .dm I is found to be equal 
to 0.69 ± 0.09 (0.76 ±0.08). 

C. Discussion 

Significant non-zero strong interaction phases are ev- 
idence that the naive factorization model is inadequate. 
Therefore, when computing the decay amplitudes, in- 
stead of using the parametrization with a\ and a 2 , the al- 
ternative parametrization in terms of isospin amplitudes 
may be more appropriate. Moreover, if we analyze the B 
decays to D^tt final states without FSI [10], we compute 
a value \a 2 \ = 0.57 ±0.07 (0.56 ±0.08). These values are 
quite different from \a 2 \ — 0.2 to 0.3 from charmonium 
final states and indicate as well the necessity of includ- 
ing strong non factorizable and process-dependent FSI 
effects in the description of B° — ► D^°h° modes. 

Various theoretical approaches that relax the condi- 
tions of naive factorization are being pursued in an effort 
to understand the emerging pattern of color-suppressed 
decay rates [10, 16, 29, 31]. 

X. SUMMARY 

We present measurements of the branching fractions 
for the color-suppressed decays B° — ► Z?°7r°, D*°7r°, D°rj, 



D*°i], D°u>, D*°lu, and D°rf . Our results are in agree- 
ment with previous measurements [2, 3] but are more 
precise. Branching fractions for B° — ► D*°rj, D*°lu, and 
D°rj' are measured for the first time with more than five- 
sigma statistical significance. We also set an upper limit 
on the branching fraction for the D*°rf decay. 

All measured color-suppressed decays have simi- 
lar branching fractions with central values between 
1.7 x 10~ 4 and 4.2 x 10~ 4 . They are all significantly 
larger than theoretical expectations based on naive 
factorization and therefore present a challenge for the 
theoretical interpretation. These results strongly suggest 
the presence of final-state re-scattering effects. 
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